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■ ABSTRACT 

O 

o . 

CN . We estimate the highest energy of proton diffusively accelerated by shock in 

W), knot Al of the jet in luminous nearby quasar 3C 273. Referring to the recent 

^ . polarization measurements using very long baseline interferometry (VLBI), we 
consider the shock propagation across magnetic field lines, namely, configuration 



X 



of the oblique shock. For larger inclination of the field lines, the effects of 
particle reflection at the shock front are more pronounced, to significantly 



in 

> , 

^ ! increase acceleration efficiency. The quasiperpendicular shock turns out to be 

needed for safely achieving the proton acceleration to the energy above 100 EeV 
^ ! (10^° eV) in a parameter domain reflecting conceivable energy restrictions. 
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: 1. INTRODUCTION 



Core-dominated quasar 3C 273 is one of the most confirmed objects on account of high 
optical luminosity and low redshift {z = 0.158; Bridle & Perley 1984). The image obtained 
with the refurbished Hubble Space Telescope (HST) revealed that the narrow optical jet 
consists of discrete knots (Bahcall et al. 1995, hereafter B95) associated with shocks. In 
radio bands the high linear polarization and featureless spectra imply that the observed 
emission from the knots is originated in the synchrotron radiation of electrons (and possibly 
positrons). If the optical emission is also of the synchrotron origin, in situ acceleration of 
electrons must be just taking place, as they are compatible with the synchrotron lifetime 
inferred from magnetic field strength (Meisenheimer et al. 1989: M89). Polarization data 



^Department of Electrical and Information Engineering, Kinki University Technical College, Mie 519- 
4395, Japan; yasuko@ktc.ac.jp. 

^Plasma Astrophysics Laboratory, Institute for Global Science, Mie 519-5203, Japan. 



- 2 - 



in optical bands are similar to those in radio ones, and suggest that the optical emission is 

attributed to the synchrotron losses of accelerated electrons (Roser & Meisenheimer 1991: 
RM91; Conway et al. 1993: C93; Roser et al. 1996). In the brightest knot A, the observed 
power-law spectrum (Si^ oc z/~") can be fitted by the index of a = 0.45 ± 0.10 (M89), close 
to the canonical value of a = (etc ~ l)/2 ~ 0.5 for the electron spectral index cte ~ 2, 
expected for Fermi acceleration at strong nonrelativistic shocks. 

As far as such an acceleration mechanism works for electrons, the same mechanism will 
operate for acceleration of ions likewise, providing their abundance in the jet is finite (e.g., 
Rawlings & Saunders 1991). One of the most promising mechanisms is the diffusive shock 
acceleration (DSA) involving resonant scattering of injected particles, in which the particles 
trapped in a mean magnetic field are resonantly scattered by the Alfven waves superposed 
on the mean field (for review, Drury 1983; Blandford & Eichler 1987; Jones & EUison 1991; 
Longair 1992: L92), so as to migrate back and forth between the upstream and downstream 
region of the shock. As a consequence, the particles are accelerated, receiving energy from 
the shock. In this context, Biermann & Strittmatter (1987) have discussed the possibility of 
ultrahigh energy acceleration of protons. However, their description was limited to a simple 
case for parallel shock propagating along magnetic field lines. Rachen & Biermann (1993) 
estimated the highest proton energy as ~ 10^^ eV, but it seems to be rather optimistic. 

For precise evaluation of the highest energy, it is important to consider the actual 
magnetic field geometry in the acceleration site. With regard to this point, a helical pattern 
of the magnetized jet has recently been observed by the Faraday rotation measures in the 
Highly Advanced Laboratory for Communications and Astronomy VLBI Space Observatory 
Programme (Asada et al. 2002), and it appears to self-organize double helical structure 
(Lobanov & Zensus 2001: LZOl). In such a filamentary jet, expected is that huge current 
(~ 10^^ A: C93) generates substantial magnetic fields transverse to the filaments (Honda, 
Meyer-ter-Vehn & Pukhov 2000; Honda & Honda 2002: HH02). Anyhow, when allowing 
shock propagation along the jet concomitant with non-parallel components of the magnetic 
fields, the shock likely crosses the field lines. This corresponds to the configuration of 
the "oblique shocf^\ for which acceleration efficiency is higher than that for the parallel 
shock, owing to the effects of mirror-refiection of particles at the shock front (Jokipii 1987; 
Ostrowski 1988: 088). 

In this Letter, we report a result of the numerical analysis for the highest energy of 
proton diffusively accelerated by the oblique shock in knot Al of 3C 273 jet. The present 
model of the oblique DSA is on the basis of the first-order Fermi mechanism including 
the resonant scattering of particles by the Kolmogorov turbulence. Concerning quench of 
acceleration and energy loss processes, more sohd argument is expanded. We demonstrate 



-3- 



that for an expected range of radiation energy density, the maximum proton energy 
exceeding 10^° eV is safely achieved, especially for quasiperpendicular shock acceleration. 

2. THE MODEL OF DIFFUSIVE SHOCK ACCELERATION AT 3C 273/ Al 

2.1. Physical Properties of Knot Al 

Below we outline the physical properties and parameters at the knot in the leading edge 
of the jet of 3C 273, which is specified as "knot Al" (B95), or "knot A" in earlier literatures. 
The knot Al is well-resolved in optical wavelengths (0.1" for HST), and its radius is about 
1 kpc (Roser et al. 2000: ROO). This optical core is encased in the larger radio cocoon 
whose half-width is ~ 2 kpc (B95); and the overall structure seems to correspond to the 
"hot spot" (in the nomenclature for FanarofT- Riley II sources; Rachen & Biermann 1993), 
which may be related to the upper limit of the size of shocks. The key feature of this 
knot is that X-ray flux is more prominent than that observed in the other knots (Marshall 
et al. 2001: MOl). Extrapolating the radio-to-optical spectrum to high-energy region 
approximately reproduces the observed X-ray continuum (ROO). Accordingly, for the X-ray 
as well, the power-law spectrum (with index a — 0.60 ± 0.05; MOl) can be explained by the 
electron-synchrotron model (Sambruna et al. 2001), rather than the other models invoking 
the inverse Compton scattering of electrons. The size of the X-ray halo, where a large 
amount of energetic particles is probably drifting, appears to larger than that of the radio 
hot spot (ROO; MOl), implying that the spatioscale of confinement region of accelerated 
particles tends to be larger than that of shock accelerator. 

In the knot distant more than ~ 20 kpc from core (B95), the flow is decelerated to the 
weak rclativistic speed of (0.21 ± 0.04)c (M89), where c is the speed of light, although in 
the vicinity of core, superluminal motion of ejecta is observed. The flow may be coUimated 
by magnetic fields (HH02), resulting in no significant radial expansion. Polarization 
measurements suggest the ordered magnetic field, and its strength is of the order of 
sub-milligauss at this knot (M89; ROO). The degree of polarization increases toward the 
direction of the bending extended structure, called "inner extension" (RM91), indicating 
that the field lines are inclined with respect to the jet axis. 

2.2. Particle Acceleration by an Oblique Shock 

For knot Al, presuming the shock speed to be nonrelativistic is adequate for modeling 
the DSA based on the Fermi-I mechanism (Gaisser 1990: G90). We consider the fast-mode 
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oblique shock, for which magnetic field strength is boosted to B2 — y'cos^'^H-^^^'sin^^Si, 
where the subscripts i ='1' and '2' indicate the upstream and downstream region, 9i denotes 
the inchnation angle of mean magnetic field line with respect to the direction normal to 
shock surface, and r is the shock compression ratio. In contrast with the case of the parallel 
shock with 61 = 62 = 0° in which particles are deflected due solely to scattering by magnetic 
fleld fluctuations, in the 61 7^ 0° case a fraction of the particles is directly reflected at 
the shock front by the boosted field in the region '2', conforming to the conservation of 
magnetic moment. Expected here is that effect of this mirror-refiection leads to significant 
reduction of acceleration time, viz., increase of acceleration efficiency (e.g.. Kirk & Heavens 
1989: KH89). Indeed, this effect has been confirmed by the Monte Carlo simulations for a 
nonrelativistic shock (Naito & Takahara 1995; Ellison, Baring & Jones 1995). 

For calculation of energy gain of particles, we transfer from plasma rest frames of the 
regions '1' and '2' to a proper frame where electric field vanishes (de Hoffmann & Teller 
1950: HT50); then, transform physical variables back to those in the original frame (088). 
The mean acceleration time can be defined as the cycle time for one back-and-forth between 
the region '1' and '2' divided by the energy gain per encounter with the shock (G90). An 
improved calculation including the mirror-reflection effects yields the following resultant of 
the mean acceleration time (Kobayakawa, Honda & Samura 2002: KHS02): 

Sr,,apac rrjg \ , sin^^i , r cos^ + p sin^ gi/(l + 77^] | 



Ul r-ll 1 + 772 (cos^^i + r^sin^^OVs 

for arbitrary species of particle 'a'. Here, t/i is the shock speed, = Va/c ~ 1, and Va the 
speed of the particle, and furthermore, 770 = ('a,\\/fg,a-i where la,\\ the mean free-path (m.f.p.) 
along the magnetic field line, r^^a — PalafnaC^ / {\(la\B\) the gyroradius in the region '1', 
7a = (1 — Pl)~^^^, and rria and Qa the rest mass and charge of the particle, respectively. Note 
the allowable range of the fleld inclination angle of 61 < ^i.max = cos^^([/i/c), where ^i,max 
is called the de Hoffmann-Teller (HT) hmit for oblique shocks (HT50; KH89; KHS02). 

In equation (1), the shock compression ratio is fixed to r = 4 for a strong nonrelativistic 
shock. We compare the magnetic fiuctuations involved in the shock to the Kolmogorov 
turbulence establishing the spectral intensity of I{k) oc k~^^^, where k is the wavenumber 
of the Alfven waves (e.g., Biermann & Strittmatter 1987: BS87). In the circumstances that 
when the wavc-particle resonance condition, Vg^a ~ k~^, is satisfied, the pitch-angle scattering 
of particles becomes effective, the m.f.p. can be denoted as ia,\\ ~ [3rg^a/(26)](rg^inax/''"g,a)^^^- 
Here, rg^max defines the maximum resonant gyroradius, and b {< 1) the ratio of turbulent to 
mean magnetic energy density. Notice that setting to 6 ~ 0(1) makes the local inclination 
angle fiuctuate around the average 9i. Nevertheless, this choice is feasible, as long as the 
trajectory of the guiding center drift of the particles bounded by the mean field is not 
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disturbed, that is, the characteristic timescale of the drift in the shock vicinity is shorter 
than the coherence time for the resonant scattering. As shown later, this condition turns 
out to be fairly satisfied in the allowable range of rg ^ < rg max- In this aspect, let us set to 
the critical value, 6=1, adequate for simple estimation of the achievable highest energy of 
accelerated particles, and convenient for making a direct comparison with a previous result 
in the special case for parallel shock [BS87, expecting 3b{Ui/c)'^ ~ 1]. 

As the jet seems to be nonuniformly filled, having helical filaments (LZOl), we allow 
the hot spot with its radius i?HS to contain multipartite shock disks, whose radial size each 
is associated with the order of maximum turbulent wavelength, ~ k^l^- In order for the 
resonance to be locked in phase, the gyroradius cannot exceed ~ k^liv which is limited by 
Rhs- That is, we have the relation of rg^a ^ max kmin ^ -Rhs- In a specific case that the 
knot feature is identified with a single shock, we reproduce k^l^ ~ -Rhs, 

i.e., Tg^max ~ -RhS) 

referred to as the Hillas criterion (Hillas 1984: H84). Below we regard, in equation (1), 
'"g,max as a variable in the range of ^ Rus- 

2.3. The Energy Constraints 

By equating the mean acceleration time (1) with the shortest timescale for the most 
severe energy restriction, we can derive the maximum possible energy of particles, defined 
as £^a,max = 7a'^aC^- For protou (a='p'), the time-balance equation can be expressed as 

^p,acc(Tp) ^1) '"gjinax) — niin[tp^syii('yp), ^p7('yp, l^rad)) ^p,esc(Tp) '"g,niax)) ^sh] ■ (2) 

Here, tp,syn — GnnipC/ (aT'm'^jpBf) defines the cooling time for synchrotron radiation 
by the accelerated protons (L92), where (Tt is the Thomson cross section. The loss 
timescale may be rewritten as tp^syn ~ 1 x 10^^(10-'-^/7p)(0.7 mG/5i)^ s. Also, the 
inelastic collision with photons involving photopionization can be a competitive energy 
loss process. For the present purpose, we employ an approximate expression of the 
loss timescale, tpj ~ (tp^syn/200)[i?^/(87rMrad)] (for 3C 273; BS87), where Urad stands 
for average energy density of target radiation fields. The above expression scales as 
tp^ ~ 1 X 10"'^^(10"'^^/7p)(10^® erg cm^'^/wrad) s- In addition to the loss due to such 
elementary processes, the particle escape itself quenches the acceleration. The escape time 
can be estimated as tp^esc — 1-5-Rxh/[c^p,||(7p5 ''"g,max)] (II84), where i?xH -Rhs) represents 
the radius of the X-ray halo confining energetic particles (see Sec. 2.1). The timescale can 
be expressed as tp^esc ~ 8 x 10-'^^(i?xH/-^p,||)(-RxH/5 kpc) s. Furthermore, the propagation 
time of shock through the jet possibly limits the acceleration, while the radial adiabatic 
expansion is less effective for the self-coUimating jet (HH02), in contrast with spherically 
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expanding supernova remnant shocks (KHS02). The shock propagation time tgh may be 
interpreted as the age of the knot, which is crudely estimated as ~ L/Upj-op, where L 
represents a distance from core to the knot and f/prop(~ Ui) the speed of proper motion of 
the knot; to give the scahng of t^h ^ 8 x 10^^(^/20 kpc)(0.25c/[/i) s. 

3. NUMERICAL SOLUTIONS: THE ACHIEVABLE HIGHEST ENERGY 

OF ACCELERATED PROTONS 

Now, given 9i, r^^max, and Uj-g^, we self-consistently solve equation (2) for 7p. Along 
the explanations mentioned above, at the moment we choose the physical parameters of 
3C 273/Al as C/i = 0.25c (M89), Bi = 700 /xG (ROO), Rhs = 2 kpc (B95), Rxh = 5 kpc 
(MOl), and L = 20 kpc (B95). It is noted that the maximum inclination angle of magnetic 
field lines reaches ^i,max = cos~^(t/i/c) = 75.5°, and ^^2,max = tan~^(r tan^i^max) — 86.3° for 
r = 4, appropriate for referring to as "quasiperpendicular shock" . 

In Figure 1 for 9i = ^i,max (left panel) and 0° (right panel), we show £^p,max = Ip^pC^ as 
a function of rg^max (~ ^min) ) for some given values of ^^rad as a parameter. In a wide range 
of smaller rg^max and larger u^ad, -E'p,max is determined by the balance of tp,acc = tpj oc 
in equation (2) (dotted curves/lines). For much smaller u^^d, -^p,max saturates, being 
determined by tp^acc = ^p.csc oc r~^^^ in a large-rg^max region (soild curve/line), whereas in a 
small-rg,max region, by tp,acc = ^p,syn for u^ad ^ 10~^° erg cm~^ (dot-dashed curve/line). Even 
for the chosen (rather small) value of L, tsh cannot be the shortest timescale in right-hand 
side of equation (2). In the both panels, bottom side of short-dashed line: Tg^p < Tg^max 
and left side of long-dashed hne: rg^^ax ^ ^hs — 2 kpc indicate the allowable and plausible 
region, respectively. Note that a point of intersection of these two lines, which represents 
the Hillas criterion, rg^p = rg^max = -RhSi is in a marginal region above the solid curve/line. 
The region of rg^p ^ ^g^max far above the short-dashed line, i.e., rj^ 1 for 6=1, can be 
compared to the "diffusive limit" (088) which violates equation (1). 

In the case of = 0°, the values of £^p,max monotonically decrease 
to exhibit the scaling of i?p,max oc r~K^^ for the cases of tp^acc = ^P7 and tp^acc = tp.syn, 
and -Ep^max oc r"^^^^ for tp,acc = ^p.esc- These tendencies reflect the resonant scattering 
theory in which longer coherence time (i.e., longer gyroperiod) leads to lower acceleration 
efficiency. For 7^ 0°, however, £'p,max is found to be significantly enhanced. As seen in 
the case of 9i — ^i,max, the resulting boost of £'p,max is prominent in the region below the 
short-dashed line. This property suggests a picture that larger value of Tjp oc (rg^max/'^g.p)^'^^ 
in equation (1) effectively represents larger anisotropy of diffusion coefficients, i.e., for 
77p » 1 (corresponding to the "free crossing limit"; 088), » k,±, where the subscripts 
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refer to the direction of the magnetic field fine; they largely assist the mirror-refiection of 
particles (see also Figure 2 below). 

We are particularly concerned with the actual maximum energy of accelerated proton 
exceeding 10^° eV. In the left panel, hatched domain indicates the window where -Bp.max 
exceeds 10^° eV for the condition of Tg < Tg^max ^ 2 kpc. In this domain, it is found that 
-E'p,max depends largely on the value of u^a,d, which involves a large uncertainty responsible 
for observation. At this juncture, we infer the value of -Uj-ad from comparing the observed 
highest frequency of radiation with the synchrotron cut-off frequency deduced from the 
calculated maximum energy of electron. This is feasible, because the X-ray from knot Al 
is arguably ascribed to the simple synchrotron emission due to accelerated electrons 
(Sec. 2.1). For electron (a='e') the method to obtain the maximum possible energy is so 
analogous to that explained in Sec. 2.3: the time-balance equation can be expressed as 
) = min[te,syn(7e),^ic(7e,iirad)]- Here, te,syn dcuotcs the famihar electron 
synchrotron timescale, and Uc — Srricc/ (4(TKN7eMrad) the timescale for the inverse Compton 
scattering (BS87; L92), where ctkn the Klein-Nishina cross section. For 9i = ^i,max, we 
have numerically calculated at i?c,max = 7c"^cC^! and the synchrotron cut-off frequency 
by using the expression of ~ {3/16){e/'mlc')E'^^^^^Bi (BS87). As a consequence, 
consistency of the calculated with the observed result fc,obs ^ 10^^ Hz (ROO) is found to 
require u^ad ^ 10~^ erg cm~^ in the region of rg^max ^ 2 kpc (not shown in figure). Note 
that the upper limit of u^ad is close to the value inferred from the energy equipartition, 
Urad = -Bf/(87r) ~ 2 X 10-^ erg cm-=^. On the other hand, for 9i = 0° we get ;^ 10^^ Hz 
(BS87) in the same region of rg^max, such that Uc <^ i^c,obs- 

For the expected value of li^ad ~ 10~^ erg cm~^, we see in Figure 1 that at rg^max = 2 kpc, 
the actual maximum energies of proton are -Bp^max ~ 2 x 10^° eV and 1 x 10^^ eV for 
^1 = ^i,max and 0°, respectively. For the 9i — ^i,max case, -Bp^max always exceeds the 
threshold of 10^° eV in the hatched domain, taking the value of about 2 x lO^'^ eV, whereas 
for ^1 = 0°, -Ep^niax cannot exceed the threshold, even though taking the peak value of 
about 6 X 10^^ cV at rg max ~ 90 pc. It is also remarked that in a possible range of 
«rad ^ 10"^ erg cm"3, -Ep^^a^ 5 x 10^° eV can be achieved for = 6'i,ma!c- 

In Figure 2 for u^ad — 10~^ erg cm~^, the ^i-dependence of -E'p,max is shown for some 
given Tg^max as a parameter. For rg max = 0.5 kpc and 1 kpc, in the whole range of 9i, £'p,max 
are determined by the balance of tp,acc = tp^, to take their maximum values at the HT 
limit, 9i = 9i „^ax- For a larger rg max = 2 kpc and smaller 0.2 kpc, this time-balance also 
governs in a major range of 9i, except for the regions of 9i <^ 8° and ^ 49°, where ii^p,max 
are determined by tp,acc = ^esc and limited by rg,p < rg max giving a constant maximum, 
respectively. Evidently, we find that the values of -E'p,max, which are, at 9i — 0°, below the 
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threshold of 10^° eV, exceed the threshold in large-^i region. For the larger Tg^max, Ep^max 
appears to be more enhanced for variation of 9i from 0° to ^i,max- 

4. CONCLUSIONS 

Referring to the VLBI observations of helical filaments in the jet of 3C 273, we have 
estimated the maximum possible energy of protons diffusively accelerated by an oblique 
shock at knot Al. The hot spot-hke feature is regarded as multipartite shock disks, or 
to contain a single shock whose radial size can be smaller than the hot spot radius. The 
complementary calculation of maximum electron energy (£'e,max/"^eC^ ^ 10^) suggests that 
the upper limit of radiation energy density is of order u^ad ~ 10~^ erg cm~^. We conclude 
that for a possible Wrad-range, protons can be accelerated beyond 10^° eV, safely by the 
quasiperpendicular shock. The present method might be applicable for solving the problem 
of in situ acceleration of particles in the other objects, although this work ignores the effects 
of particle transport, which may be important for making a comparison with results of 
precedent (Takeda et al. 2003; Abbasi et al. 2004) and future experiments (at the Pierre 
Auger Observatory, the Extreme Universe Space Observatory, and so forth). 
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Fig. 1. — The maximum possible energy of accelerated proton E'p^max vs. the maximum 
resonant gyroradius rg^max for the inclination angles of magnetic field lines, 9i — 6'i,max = 
75.5° (left) and 9i — 0° (right). The figures have the same axes. The solutions of 
equation (2) are plotted for the balance of the mean acceleration time, ip,acc, with the loss 
timescales for escape (tp,esc; solid curve/line), synchrotron {tp^syn, dot-dashed curve/line), and 
photopionization (tp^; dotted curves/lines). For ip^acc = tp^, the curves/lines are displayed 
each factor 10 for radiation energy density, "Uradj in the unit of erg cm~^. Note that those 
plots for -Urad ^ 10^^° erg cm^^ overlap the dot-dashed curve/line. In the left panel, hatched 
domain indicates the window where £^p,max exceeds the threshold value of 10^° eV (triple 
dot-dashed line) for the conditions of Tg^p < rg^max (bottom side of short-dashed hue) and 
^g,max ^ 2 kpc (left side of long-dashed line) . 
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Magnetic field inclination (deg) 

Fig. 2. — The dependence of the maximum proton energy -Bp^max on the inchnation angle 6i 
for the parameters of rg^max = 0.2 (dotted), 0.5 (dot-dashed), 1 (dashed), and 2 kpc (sohd). 
Here, we have set to Wrad = 10~^ erg cm~^. In the inset, the inchnation angles 6i are limited 
by the maximum values of ^i^max = 75.5° for the shock speed of 0.25c, and 6'2,max = 86.3° for 
the shock compression ratio 4. 



